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Black holes generate coUimated, relativistic jets which have been observed in 
gamma-ray bursts (GRBs), microquasars, and at the center of some galaxies 
(active galactic nuclei; AGN). How jet physics scales from stellar black holes in 
GRBs to the supermassive ones in AGNs is still unknown. Here we show that 
jets produced by AGNs and GRBs exhibit the same correlation between the 
kinetic power carried by accelerated particles and the 7-ray luminosity, with 
AGNs and GRBs lying at the low and high-luminosity ends, respectively, of the 
correlation. This result implies that the efficiency of energy dissipation in jets 
produced in black hole systems is similar over 10 orders of magnitude in jet 
power, establishing a physical analogy between AGN and GRBs. 

Relativistic jets are ubiquitous in the cosmos and have been observed in a diverse range of 
black hole systems spanning from stellar mass (~ IOMq) to supermassive scales (~ 10^ — 
10^°Mq), in particular in the bright flashes of gamma-rays known as GRBs (i, 2), the miniature 
versions of quasars lurking in our galaxy known as "microquasars" (5) and AGNs (4,5). Despite 
decades of observations at almost all wavelengths and considerable theoretical efforts, there are 
still many aspects of black hole jets which remain mysterious: the mechanism(s) responsible for 
their formation and the nature of their energetics as well as their high-energy radiation (6, 7). Jets 
and outflows from supermassive black holes have important feedback effects on scales ranging 
from their host galaxies to groups and clusters of galaxies (8). Hence, a better understanding 
of the physics of jets is required in order to have a more complete picture of the formation 
and evolution of large-scale structures in the universe and the coevolution of black holes and 
galaxies (9). 

One outstanding question is how the jet physics scale with mass from stellar to supermassive 
black holes. Interestingly, there is evidence that jets behave in similar ways in microquasars and 
radio-loud AGN (10-12). However, a clear connection between AGN and GRBs has not been 
established yet, although recent work provides encouraging results (13, 14). 

As a first step in understanding how the properties of jets vary across the mass scale, we 
focus on the energetics of jets produced in AGNs and GRBs. Therefore, we searched the litera- 
ture for published and archival observations that allow us to estimate the jet radiative output and 
the kinetic power for a sample of black hole systems in which the jet is closely aligned with our 
line of sight and characterized by a broad range of masses. For this reason, our sample consists 
of blazars - AGNs with their jets oriented toward Earth (15) - and GRBs, the spectral energy 
distributions of which are completely dominated by the jet due to beaming effects. 

We used as a proxy of the jet bolometric luminosity the observed 7-ray luminosity 
which is isotropically equivalent. In order to estimate the kinetic power Pjet, we use extended 
radio luminosities for the blazars whereas for the GRBs we relied on the afterglow measure- 
ments in radio or X-rays. Therefore, the availability of these observables restricted our sample 
to 234 blazars (106 BL Lacs and 128 flat-spectrum radio quasars - FSRQs; see Table SI) and 
54 GRBs (49 long and 5 short GRBs, all with known redshifts z; see Table S2). For blazars, 
was estimated from the 7-ray energy flux and the spectral index measured with Fermi Large 
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Area Telescope (LAT) (16); Pjct was estimated using an empirical correlation which relates the 
Very Large Array (VLA) extended radio emission and the jet kinetic power (17, 18). For GRBs, 
2^ ISO _ ^iso(^]^ _j_ z)/t<^Q where tgo is the burst duration and E^^° is the isotropically equivalent 
energy radiated during the prompt emission phase and measured with different telescopes (21 
observed with either BeppoS AX, BATSE, HETE, HETE-2 or Integral, 24 with Swift Burst Alert 
Telescope - BAT - and 10 with Fermi). Pjet was computed as Pjet = fbE]^°{l + z)/tgQ where 
E^^° is the kinetic energy estimated from the radio (VLA) or X-ray (Chandra) luminosity during 
the afterglow phase using the standard afterglow model (19), = 1 — cos 6 is the "beaming 
factor" and 6 is the radiation cone half-opening angle which is the same as the jet opening angle 
estimated from the jet break in the GRB afterglow lightcurve (20). 

We first compared the relative trends of and Pjct for the blazar and GRB population 
separately (Figure 1). The Pearson correlation coefficients of 0.85 and 0.8 obtained for blazars 
and GRBs respectively, indicate a strong correlation within each group of sources. However, 
the U^°-P-^ct trend is different for GRBs and blazars as shown by the fits to the data (Fig. 1). 

We computed the intrinsic luminosity, L, for GRBs and blazars by correcting for the 
opening angle or beaming factor, /b, such that L = For GRBs, the beaming factor 

is computed from the jet opening angle 9j as 1 — cos6'j (22); for blazars, is estimated as 
1 — cos 1 /r where T is the bulk Lorentz factor of the flow, since AGNs obey 9^ < 1/T (23, 24). 
While an estimate of Oj is available for each GRB in the sample, T is only available for a subset 
of 41 blazars. Figure 2 shows an anti correlation between and fb for both GRBs and blazars 
with compatible indices when fit with a power law. Because 9 is not available for the whole 
blazar sample, we used the power-law fit of vs as an estimator for fb. 

As with U^" and Pjct, L and Pjct are strongly correlated within the GRB and AGN samples 
(Fig. 3). However, they follow the same trend within the narrow uncertainties and the whole 
GRB and blazar sample can be fit adequately with a power law over 10 orders of magnitude 
in luminosity. Therefore, the relativistic jets in GRBs and blazars are consistent with obeying 
the relation Pjet ~ 4.6 x 10^^ (L/10^^)°'^* erg s~^, within the measurement uncertainties. In 
other words, once "black hole engines" produce relativistic jets, they seem to do so maintaining 
the same coupling between the total power carried by the jet and power radiated away. This 
universal scaling for the energetics of jets is maintained across the mass scale regardless of the 
different environments and accretion flow conditions around the compact object. 

Figure 4 indicates that most of the jets in our sample dissipate at least 3% of the power 
carried by the jet as radiation and overall they can radiate as much 15%. This range of ef- 
ficiencies is considerably higher than previous estimates for AGNs based on radio to X-rays 
luminosities (25, 26) but they are in agreement with results obtained from blazar broadband 
spectral models (27, 28) as well as GRB afterglow studies (29-31). Efficient heating of elec- 
trons seems to be a universal property of relativistic magnetized shocks according to numerical 
simulations (32) which demonstrate that electrons retain > 15% of the pre-shock energy. If 
most of the post-shock energy is radiated away, these theoretical results could pave the way to 
an understanding of the high dissipation efficiencies that we find. 

Our results suggest that there is a single fundamental mechanism to produce relativistic jets 
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in the Universe. The analogy known to exist between microquasars and AGNs (3, 10, 11) can 
be extended to the gamma-ray bursts with the fundamental difference that whereas AGNs and 
microquasars undergo recurrent activity, GRBs experience only one episode of hyperaccretion. 
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Figure 1: The relation between the jet kinetic power and the isotropically-equivalent 7-ray 
luminosity for AGNs and GRBs. Error bars, la. We fitted the two populations separately 
using a symmetric least-squares method (orthogonal BCES with bootstrapping; 21). The blazar 
and ORB best-fit models corresponds to the solid and dashed lines, respectively (logPjot = 
A\ogU^° + B). The best-fit parameters obtained for the blazars are A = 0.51 ± 0.02 and 
B = 21.2 ± 1.1; for the GRBs, A = 0.74 ± 0.08 and B = 11.8 ± 4.1. The scatter about the 
best-fit is 0.5 dex and 0.8 dex for the blazars and GRBs, respectively. The 2a confidence band of 
the fits is shown as the gray shaded regions (barely visible for blazars). The two correlations do 
not agree at > 5a level. We also include for illustration XRF 020903 and GRB 090423 (yellow 
circles) as well as the two recent tidal disruption flares (TDFs) detected with Swift which are 
presumably due to the onset of relativistic jets from the tidal disruption of stars by supermassive 
black holes (36). We do not consider these sources in the statistics since we only have limits on 
their luminosities 
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Figure 2: The relation between the apparent 7-ray luminosity and the beaming factor for blazars 
(left panel) and GRBs (right panel). We find r = —0.53 and —0.56 for blazars and GRBs, 
respectively, indicating anti correlations significant at the 3.6a and 4.4cr levels respectively. The 
solid lines correspond to the best-fit linear models obtained with the symmetric least-squares fit 
and are given by ^ 5 x 10-^(L^9°)-°-39±o i5 and ^ 0.03(Li9°)-°-24±o o6 for blazars and GRBs, 
respectively. The gray shaded region corresponds to the la confidence band and the blue and 
yellow regions are the la prediction bands, which quantify the scatter about the best-fits. 
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Figure 3: The relation between the collimation-corrected 7-ray luminosity L = fbL^^° and the 
kinetic power for AGNs and GRBs. The shaded regions display the 2a confidence band of 
the fits. The blazar and ORB best-fit models (dashed and dotted lines, respectively) follow 
correlations which are consistent, within the uncertainties, with the best-fit model obtained 
from the joint data set (solid line). In other words, using L instead of U^'^ leads to correlations 
for AGNs and GRBs which are consistent with each other (compare to Fig. 1). The best-fit 
parameters obtained from the combined data set are a = 0.98 ± 0.02 and /3 = 1.6 ± 0.9 
where logPjet = alogL + (3. The scatter about the best-fit is 0.64 dex. The yellow data 
points correspond to XRF 020903 and GRB 090423, which we do not take into account in the 
statistics. 
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Figure 4: Distribution of the Icr lower limits on the jet radiative efficiency (the fraction of the 
total jet power which is converted to 7-rays) erad = L/{L + Pjct) for AGNs and GRBs. The 
vertical solid lines indicate the median values of the lower limits and the dashed lines represent 
the median values of e^ad for each sample. Most of the sources are characterized by Crad > 3%. 
The median efficiencies correspond to about 15%, keeping in mind that these estimates are 
affected by ~ 0.5 — 0.7 dex uncertainties on average. 
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Supplementary Materials 



Data 

Table SI lists the properties of the 234 blazars in our sample (106 BL Lac objects and 128 
FSRQs). Our sample includes the well-studied blazars OJ 287, 3C 454.3, 3C 279 and 3C 273. 
The redshifts were estimated as described in (17). 

Table S2 lists the corresponding properties of the 54 GRBs (49 long and 5 short). We 
include in our data the sub-energetic bursts GRB 031203 (57) and GRB 980425 (38). The latter 
one is associated with the nearby (distance ~ 40 Mpc) supernova 1998bw (39, 40). We also 
include in our sample the naked-eye GRB 0803 19B (42). We include in Table S2 the X-ray 
flash (XRF) 020903 (41) and the most distant cosmic explosion ever detected, GRB 090423 at 
z a; 8.2 (43, 44) but we do not take these GRBs into account in the statistics, since we have only 
limits on their coUimation-corrected energetics. 



Calculation of 7-ray luminosity and jet power 

Blazars 

7-ray luminosity: 

In order to calculate the total Fermi 7-ray luminosity, we follow a procedure similar to that 
of (46) (their equation 1). The procedure to calculate the /c-corrected band luminosity depends 
on the type of model used in the 2FGL analysis (16, 47). In the case of a power law energy 
model, the total 100 MeV to 100 GeV A; -corrected luminosity is calculated from the energy flux 
(S^) given in the catalog: 

I}-- ^ Ai^dl—^^ , (1) 

where c^l is the luminosity distance in cm^ and is the (energy) spectral slope over the whole 
band. 

In the 2FGL catalog, some blazars are now modeled with the 'Log Parabolic' form. The 
/c-corrected energy flux in this case must be calculated numerically. The integral form is 

where -Ei = 0.1 GeV and E2 = 100 GeV and we have used the fit given in the 2FGL catalog for 
each source, with values a. (column name 'spectralJndex'), (3 ('beta'), and Eq ('pivot_energy'), 
and K ('flux.density'). The constant x = 1.6 erg MeV GeV"^ gives the energy flux in final units 
of erg cm~^ s~^. 

The band luminosity for the Log-Parabolic model case can then be simply calculated from 

L-° = 47r45;. (3) 
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We calculate the uncertainty in propagating the error associated with and quoted 
in the 2FGL. The average uncertainty in corresponds to 0.05 dex. 

In order to estimate the uncertainty affecting the coUimation-corrected luminosities L we 
first evaluate the error in the beaming angle Q or correspondingly F. The uncertainty in F is 
dominated by the uncertainty in the variability Doppler factor whereas the uncertainty in the 
apparent speed does not contribute significantly to the error budget of F. The uncertainty in 
the variability Doppler factor is ~ 27% (1 s.d.; 4S). Therefore, for the blazars with direct 
estimates of F available, the relative uncertainty in F is 0.3 (24, 48) which translates to an 
average uncertainty of 0.26 dex in L for these blazars. For the blazars without direct estimates 
of F, we estimate the uncertainty in L using the prediction band of the — fb relation shown 
in Fig. 2. The plotted prediction band corresponds to a relative uncertainty of 0.69 in 9. The 
resulting average uncertainty affecting L for the blazars without direct estimates of F is then 0.6 
dex. 

Kinetic power: 

Following (77), we estimate the jet kinetic power by using the correlation between the ex- 
tended radio emission and the jet power (18, 49). Cavagnolo et al. searched for X-ray cavities in 
different systems including giant elliptical galaxies and cD galaxies and estimated the jet power 
required to inflate these cavities or bubbles, obtaining the tight correlation 



between the "cavity" power and the radio luminosity. Hence, assuming Pjct = -Pcav we can 
estimate the jet kinetic power for the blazars which have extended radio emission observed with 
theVLAa?). 

The uncertainty in Pjct is dominated by the scatter in the correlation of (18) and corresponds 
to 0.7 dex. 



^-ray luminosity: 

The current scenario for GRBs (50, 51) posits that initially most of the energy produced by 
the GRB is in kinetic form produced during the short "active" state of the stellar-mass central 
engine. A certain fraction of initial energy is converted after a few seconds mostly to 7-rays 
observed during the prompt emission, by means of internal shocks in the jet (52). The ultrarel- 
ativistic jet produced in the explosion later on collides with the circumburst medium producing 
the afterglow. Two crucial quantities which we use in this work are the radiative and kinetic 
energies released by the GRBs during their short period of activity. 

The isotropically equivalent energy radiated in 7-rays E]^° is directly available from mea- 
surements. It was measured for the GRBs using a variety of different telescopes including 
pre-Swift telescopes (BeppoSAX, BATSE, HETE, HETE-2 and Integral) as well as Swift and 




(4) 



GRBs 
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Fermi (see Table S2). We calculate the isotropically-equivalent 7-ray luminosity as 



^90 

where tgo is the duration containing 90% of the fluence in the observer frame. 

The energy range in which the fluence is measured is typically ~ 10 keV - 10 MeV. Most 
of the radiative energy released in the GRB jet during the prompt emission is contained in this 
energy range according to the latest GRB SEDs observed (55). 

We adopt an uncertainty of 0.2 dex on the values of which corresponds to the typical un- 
certainty affecting the GRBs in the sample studied by {54). Therefore, the resulting uncertainty 
in the value of L'*^" corresponds to 0.2 dex. 

The coUimation-corrected 7-ray luminosity is computed as 

L = ^^^A^;^ (6) 

190 

where is the beaming factor (/^ = 1 — cos ^j) and is the jet half-opening angle. This relies 
on the afterglow lightcurve displaying a jet break which is used to estimate Q (22). (9^) ^ 8° 
for the sample and (/(,) ^ 9 x 10~^. We adopt an uncertainty of 0.1 dex on the values of 6'jet, 
which corresponds to the typical uncertainty in the values of 9j for the sample studied by (54). 

We calculated the uncertainty in L using error propagation from the uncertainties in 
and 9, obtaining that the uncertainty in L is fa 0.3 dex. 

Kinetic energy: 

The jet kinetic energy is estimated from the radio or X-ray afterglow lightcurve using the 
fireball model (19). For most of the GRBs, X-ray data were used to determine this energy 
(31, 54). In a few cases, radio data were used (35, 37, 41). 

The standard fireball afterglow model depends on five model parameters: the explosion 
kinetic energy E^,^°, the density of the circumburst environment n (with which the jet collides), 
the spectral index of the electron energy distribution p and the fractions eg and of shock 
thermal energy carried by electrons and magnetic field, respectively. The typical values adopted 
for these parameters are p = 2.2, n — 1 cm~^, €e — 0.1 and €b — 0.01 (e.g., 19, 31). The 
afterglow model relates the specific flux at a certain frequency and at a specific time after the 
burst (typically 10 hours in the observed frame) with the kinetic energy. 

The measurement of E]^° can be impacted by the afterglow plateau which possibly corre- 
sponds to a late activity of the central engine. Indeed, (30) demonstrated the impact of choosing 
two different times for the measurement of E^^°: tdec (deceleration time) or tb (injection break 
time). Adopting tdec would lead to an underestimation of E^^° because this choice of time does 
not include the plateau. 

The kinetic energy estimates we used in our analysis were computed at either ~ 10 h or 
~ 24 h (cf. Table S2). We verified that these times are usually beyond the "break time" of the 
plateau reported in (45). Therefore, the measurements of E^^° used in this paper correspond to 
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conservative estimates. Moreover, the choice of the time at either ~ 10 h and ~ 24 h affects 
very little the measurement of i?^^" (31). 

The typical uncertainty in E^^" due to observational errors is ~ 0.3 dex (30, 54). The value 
of E^^° is also sensitive to the values of parameters which regulate the microphysics of the 
fireball afterglow model and are poorly constrained. The systematic error affecting E'^'^ due to 
the uncertainties in the parameters eg and eg of the fireball afterglow model can be as high as 
0.45 dex (29, 30). Therefore, we combined the uncertainty resulting from the observational and 
systematic sources of errors in quadrature and conservatively adopt an uncertainty of 0.5 dex 
for E'^"". 

Kinetic power: 

The jet kinetic power is computed as 



D^" and Pjct should be thought as the average luminosities over the duration tgo of the prompt 
emission phase, i.e. the average luminosities over the timescale during which the central engine 
is producing the jet. We calculated the uncertainty in Pjet using error propagation from the 
uncertainties in E^ and 9, obtaining that the uncertainty affecting Pjet is 0.54 dex. 

Linear regression method 

Here we present more details about the linear regression method that we use in the paper. 

We fitted the datasets using the BCES (bivariate correlated errors and intrinsic scatter) re- 
gression method (27) which takes into account measurement errors in both the "X" and "F" 
coordinates and the intrinsic scatter in the data. This method has been widely used in fitting 
datasets in the astronomical community (70, 71). 

It is not clear in the data sets analyzed in this work which quantities should be treated as 
the dependent variables and which ones should be treated as independent from a physical point 
of view. There is no a priori reason to expect the luminosity to be the independent variable 
as opposed to the kinetic power (or beaming factor). For this reason, we treat the variables 
symmetrically and adopt the BCES orthogonal regression method, which minimizes the squared 
orthogonal distances. Uncertainties on the parameters derived from the fits are estimated after 
carrying out 100000 bootstrap resamples of the data. 

Partial correlation analysis 

When studying correlations between luminosities one should be careful to take into account 
their common dependence on the distance {11). We performed a partial correlation analysis of 
the common dependence of L'^° and Pjet on the distance using the partial Kendall's r correlation 




(7) 



test (72). 
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Objects 


N 




a 


-Pimll 


Signif. rejection null 


Blazars 


234 


0.3 


0.04 


5 X 10-^5 


7.8 


GRBs 


54 


0.4 


0.08 


1.6 X 10"^ 


5.2 



Table S3 Results of partial correlation analysis with X — Pjet, Y — L'®° and Z — log (II. 
Notes. Column (1): subsample. Column (2): Number of sources. Column (3)-(6): results of 
partial correlation analysis; r is the partial Kendall's correlation coefficient; a is the square root 
of the calculated variance; Pnuu is the probability for accepting the null hypothesis that there 
is no correlation between X and Y; Column (6) gives the associated significance in standard 
deviations with which the null hypothesis is rejected. 



Objects 


N 


r 


a 




Signif. rejection null 


All objects 


288 


0.57 


0.06 


< 10-10 


9.4 


Blazars 


234 


0.24 


0.03 


2 X 10-12 


7 


GRBs 


54 


0.47 


0.08 


10-9 


6.1 



Table S4 Same as Table S3 with X = Pjet, F = L and Z = log d^. 

We applied this test to our data considering X = Pjct, Y = U^" and Z = log d^. Table S3 
lists the results of the partial correlation analysis. This test demonstrates that the p-value of the 
null hypothesis (i.e. no correlation between X and Y) is 1.2 x 10-^ when considering the GRB 
subsample and < 10~^^ when considering the blazar subsample or the combined blazar-GRB 
sample. Therefore, the — Pjct correlation is strong and not a distance-driven artifact. 

We then applied the partial correlation test to the data considering X = Pjct, Y = L and 
Z = logdi, i.e. we replaced the isotropically-equivalent luminosity in Y with the coUimation- 
corrected one. Table S4 lists the results of the analysis. This test demonstrates that the p-value 
of the null hypothesis is <^ lO-^^, fa lO-^^ and IQ-^ when considering the combined blazar- 
GRB sample, the blazars and the GRB, respectively. Hence, the L — Pjet correlation remains 
very strong after correcting for beaming. 

Blazar luminosity estimates: Impact of the synchrotron peak 

We discuss in this section the impact of the lower energy synchrotron peak - observed in the 
spectral energy distribution of blazars (74, 75) - in estimating the jet radiative luminosity. 

A subset of 13 1 blazars - roughly half of the original AGN sample - have adequate sampling 
of their multiwavelength spectral energy distributions which allow us to quantify the impact of 
the synchrotron peak in the estimate of the jet radiative luminosity. We estimated the bolometric 
luminosity as U^°i = L'^° + Lf°^ where Lf°^ is the isotropically-equivalent luminosity of the 
synchrotron peak. We obtained that the values of for the quasars are not much different 
from L"^° (on average by a factor of ^ 1.7), whereas the values of Lj^°j for the BL Lacs are 
somewhat different (on average by a factor of pa 2.8). 
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We computed the intrinsic bolometric luminosity as Lboi = /b-^^boi' correcting for the 
opening angle or beaming factor. As previously discussed, estimates of fb for blazars rely on 
measurements of T which are not available for all the blazars in our sample. For this reason, 
we use the anti correlation between and ff, as an estimator of fh for the blazars without 
measurements of T. 

We show in Figure SI the resulting relation between Lboi and Pjct compared to the L — Pjet 
fit derived before (cf. Fig. 3). Figure S 1 illustrates that the Lboi — -Pjet and L — P^^t best-fits are 
characterized by very similar slopes. However, the fit based on Lboi is characterized by slightly 
higher radiative efficiencies: for a given jet power, the jet luminosity is higher on average by a 
factor of ~ 2 compared to the fit based on L. Therefore, the results based on Lboi strengthen 
our conclusion that AGN jets have high radiative efficiencies while still being in qualitative 
agreement with the GRB result. 
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Figure SI The relation between the coUimation-corrected bolometric luminosity Lboi and the 
kinetic power for the 131 blazars with measurements of synchrotron peak luminosity. The 
solid line in this figure corresponds to the fit based on Lboi whereas the dashed line is the fit 
based on L. The shaded regions display the 2a confidence band of the fits. Error bars, la. 
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